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In our recent paper we presented empirical evidences suggesting that electrons in cuprate su- 
perconductors are strongly coupled to unscreened c-axis polar phonons. In the overdoped regime 
the c-axis metallizes and we present here simple theoretical arguments demonstrating that the ob- 
served effect of the metallic c-axis screening on the polar electron-phonon coupling is consistent 
with a strongly overdamped c-axis charge dynamics in the optimally doped system, becoming less 
dissipative in the overdoped regime. 
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I. INTRODUCTION 

As electrodynamical media the cuprate high-T c super- 
conductors are highly anomalous. In their normal state 
they behave like metals in the planar ab-directions while 
they can be regarded as dielectrics along the interpla- 
nar c-direction, at least at finite frequencies. Only in 
the regime where the superconductivity start to degrade 
because of too high doping levels, metallic 'normalcy' de- 
velops along the c-axisi 

Electrodynamical 'anomaly' raises a general question 
regarding the nature of electron-phonon coupling in the 
cuprates. Insulating cuprates are as any other oxide 
characterized by a highly polarizable lattice. The con- 
sequence is that one is dealing with strongly 'polar' 
electron-phonon (EP) couplings, associated with the long 
range nature of the Coulomb interaction between carri- 
ers and the ionic lattice 2 . In normal metals these long 
range interactions are diminished by metallic screening 
and one is left with rather weak, residual short range EP 
interactions. Although this is also the case in the pla- 
nar directions of the cuprates, due to incoherent electron 
motion perpendicular to the planes we argue that the c- 
axis phonons are largely unscreened in this direction and 
strong polar EP-couplings are active 3 . 

In our recent paper—, angle-resolved photoemission 
(ARPES) studies of optimally-doped and overdoped 
Bi2Sr2Cu06 (Bi2201) superconductors have indicated 
that polar c-axis couplings play an important role. En- 
ergy dispersion "kinks" resulting from coupling to well 
defined bosonic modes, which we argued to be phonons 
rather than spin resonance modes, are weaker at low 
binding energies in overdoped compared to optimally 
doped Bi2201. This case is twofold: (a) a good descrip- 
tion is obtained for the electronic self-energy in the nodal 
momentum direction in an optimally doped single layer 
Bi2201 superconductor by assuming that the Eliashbcrg 
function a 2 F(u>) largely coincides with the measured c- 
axis electron energy loss function; (b) The electronic self- 



energy drastically changes in a strongly overdoped sys- 
tem and we show that these changes are accounted for as- 
suming a frequency dependent screening characterized by 
a scale ui S cr,c — 60meV. Here we will present calculations 
further substantiating this claim. We will give the rea- 
sons why the c-axis optical loss function should be a good 
model for the Eliashbcrg function associated with the 
polar couplings in the optimally doped system. Further- 
more, we will show that the way the self-energy dimin- 
ishes in the overdoped regime is qualitatively- and quan- 
titatively consistent with what has been learned from op- 
tical measurements about the way metallicity develops in 
the overdoped regime. 

We remark that a great deal of focus has been placed 
on in-plane phonons such as the breathing modes show- 
ing large softenings with doping in the cuprates^. These 
modes have been shown to contribute to the electronic 
self-energy for nodal electrons^. Moreover, the breathing 
phonons provide signatures on the spectral function for 
the formation of charge ordering 7 and may contribute to 
d— wave pairing in the presence of strong correlations^. 
Rather than taking into account the interesting and rele- 
vant issue of planar phonon couplings, our focus here will 
be on the c— axis phonons in order to explore the role of 
screening. 

The c-axis is insulating, not because a lack of carri- 
ers, but instead due to the fact that charge dynamics is 
strongly overdamped and incoherent. To be discussed at 
the end, the damping rate is decreasing rapidly in the 
overdoped regime, and this is exactly what is needed to 
explain the behavior of the self-energy. Van der Marel 
and Kim 9 pointed out some time ago that the effective 
c-axis plasmon frequency w PjC = 4-7m c e 2 /m c (n c and m c 
are the c-axis effective carrier density and band mass, re- 
spectively) as determined by an optical sum rule can be 
as large as 1/4 of the planar plasma frequency w p , a & — 1 
eV— . The essence of overdamped dynamics is captured 
by the expression for the dielectric function associated 
with a plasmon with damping rate T and frequency lo v . c 
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at long wavelengths, 



{iJ 2 + iLuT) 



(1) 



It is demonstrated (see Eq. l(T7)l that in the overdamped 
case r >> ujp c the medium behaves like a dielectric at 
frequencies larger than the characteristic frequency, 



eoor 



(2) 



In the cuprates, it will be shown that in the doping 
range up to optimal doping, this frequency is smaller than 
the typical optical phonon frequencies with the effect that 
the c-axis optical conductivity <j\(uj) looks like that of a 
polar insulator, dominated by unscreened polar phonons. 
To get the impression of the magnitude of uj S cr,c, we as- 
sume a simple Drude form for the optical conductivity 
(see Eq. $18\i ). With this form, uj SC r,c is approximately 
proportional to the height of the electronic background in 
the c-axis cr^(w) and this background is found to increase 
rapidly in the overdoped regime in LSCO systemii. We 
find that this change is just of the right magnitude to 
explain the changes in the self-energy of the overdoped 
system. 

The justification for taking the c-axis EP couplings be- 
ing representative of the interactions showing up in the 
self-energy will be discussed first. The actual importance 
of polar EP-couplings in the cuprates is a priori a rather 
subtle affair. On the one hand, there are actually very 
few phonons which can be regarded as unscreened. The 
cuprates can be viewed as a stack of metallic sheets. For 
one limit, when considering phonons with q purely along 
c-axis (q = q c ), the screening frequency uj scr ^ c should be 
set by the effective c-axis plasmon frequency w PjC , such 
that all phonons with frequency u> p h,i(q = q c ) < ^scr.c 
can be regarded as screened. This Lo scr ,c will be later 
shown to be defined by Eq. (|17() . Similarly, for the other 
limit, when considering phonons with q purely in the ab- 
plane (q = q a b), ^>scr,ab should be set by the planar plas- 
mon frequency w p . a b- Since u> Pi ab is large compared to 
oj p < m& : this implies that phonons with q along the c-axis 
will be left unscreened compared to phonons with q in 
ab plane. In a more general case, at a 3D momentum 
Q = Qab + Qc, the effective 3D plasmon frequency w p ^ for 
q — > O 9,12 (also see Appendix A) is given by, 
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The polar coupling requires that the phonon frequencies 
exceed the 3D screening frequency, which should be set 
by this to 2 Hence, this implies that only phonons in a 
narrow cone around q a b — can contribute to the polar 
couplings. It is questionable whether this small q a b will 
show up in the polar couplings. Given the unknowns, this 
question is not easy to answer on theoretical grounds. 
On the other hand, leaving the overall strength aside, 
the frequency dependence of the electronic self-energy 



can be easily calculated 13,14 and in section III we will 
compare the experiment- with the calculation, assuming 
the predominant electron-phonon coupling is due to the 
unscreened c-axis phonons. 



II. REAL PART OF ELECTRON SELF-ENERGY 

In the following, the real part of self-energy, i?e(£), 
contributed from interactions which can be captured by 
the c-axis optical dielectric function, e c (w), will be calcu- 
lated where the strength of interactions as a function of 
energy is represented by the Eliashberg function, a 2 F(ui). 
Regarding the empirical result shown in the experimen- 
tal paper—, the main results in this paper are that a) the 
c-axis loss function, Im{— l/e c (w)), is a good representa- 
tive of a 2 F(us) for the optimally doped cuprate and b) 
when including screening effect, the expression of a 2 F(u>) 
becomes Eq. Ijlfijl giving the good agreement of i?e(S) of 
the overdoped cuprate. 

The electron self-energy can be given by (<5 — > + ), 

E(^,e) = -— — - J d 3 Pl J alto J deJmWip-p^Lu) 

(4) 
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where G and W are the retarded electron propagator 
and the effective interaction, respectively. Here we ne- 
glect the detailed momentum structure of the electron- 
phonon coupling, which although important for a subset 
of phonons (such as the Bi g and Bi„ c-axis phonons), 
is motivated by the interest in small q phonons provid- 
ing a more or less isotropic charge displacement along 
the c-axis. We neglect vertex corrections which may be 
important for polaron formation but are not crucial to 
develop the ideas of screening. 

When only considering the interaction of electrons to 
polar c-axis phonons and electrons to electrons, the ef- 
fective interaction (see Appendix B) can be written as 



W eff (q,u) 



(5) 



where e tot (q,u)) — too - v q P e (q 7 uj) + e ph (q,u) is the to- 
tal dielectric function in which the core polarization, the 
phononic- and electronic polarizations, in terms of the 
electronic polarizability P e (q,u)), enter additively. 

The effective interaction includes both Coulomb and 
EP interactions, and the manifestation of screening is 
that the interactions are mixed as Coulomb interactions 
modify the phonon propagator D and vice- versa. This ef- 
fective interaction can be rewritten a.si& {M q is the bare 
EP vertex neglecting any fermionic momentum depen- 
dence) 



W e ff(q,Cj) = 



£ei (q, w) 



,M 2 q D{gM 

<?ei{q,u) 



(6) 
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where e e i is the electronic polarizability Eq. (JTJ and 
v q = Aite 2 /q 2 ] this is a general expression for phonon- 
modulated uniform charge displacements. 

We now assert that the effective interaction We// 
is dominated by the coupling to the unscreened c-axis 
phonons. Next we assume that momentum dependences 
are smooth. This is definitely not problematic for the 
phonon-propagators hidden in W e t / because the impor- 
tant phonons are of the high energy optical variety having 
small dispersions. This is further helped by the fact that 
because of the metallic screening in the planar directions, 
only phonons with a small planar momentum q a b can be 
regarded as unscreened. A more delicate issue relates to 
the electron momentum dependence of the bare vertex 
having implications for how the self-energy will depend 
on p. Elsewhere we will show that this can give rise to 
highly interesting effects associated with the screening 
physics^ but here we will focus on the frequency depen- 
dences in a given momentum direction. 

Neglecting momentum dependences, the phonon po- 
larizability can be straightforwardly parameterized, by 
a sum over phonons with transversal frequencies u>Tj, 
damping -fj and strengths Sj, as tabulated by Tsvetkov 
et ali&, 

2 

^»h=E 2 Sj 7L- ( 7 ) 

Taking for the electronic polarizability v q P e (uj) = 
lo 2 ^ c /u) 2 and a plasmon frequency w PiC large compared to 
all phonon frequencies one recovers the standard results 
for the EP coupling in metals, but we now focus on the 
category of unscreened phonons, such that v q P e (uj) = 0. 

In the insulator, the phonon propagator turns into 
the bare phonon propagator D — > D° which becomes 
in terms of the parametrization Eq. JJJ 

Coo + £ph\U) 

where =v q /e OQ . 

The Eliashberg function is defined as the Fermi-surface 
average— 

a 2 F t (w = -— - / s k Z k 2 ImD(k~k',uj) 

(9) 

while the real part of the electronic self-energy, calculated 
from the Eliashberg function, is given byi^, 

feS(p,e) = d^a 2 F(co-p,e)K^,^) (10) 

where K(y,y') = - f^dxffa - y)2y'/(x 2 - y 2 ) with 
f(x) being the Fermi distribution function. 



A. a 2 F(uj) for c-axis insulating system 
(optimally-doped) 

Assuming smooth momentum dependences, and ob- 
serving that in the insulator e e i = e^, one immediately 
infers that the Eliashberg function for the polar phonons 
can be written as, 

\M 2 I / 1 \ 
a 2 F(tu)=b l q ^ 0l ImD (ij) = b-Iml — (11) 

v q ~o V e to*M / 

where lumping together in the factor b is the numeri- 
cal factors coming from the momentum integrations (see 
Appendix C). Since momentum dependences are weak we 
can as well take the optically measured c-axis loss func- 
tion im(— l/e c (w)) where e c (w) = e to t(q a b = 0, q c — > 0, u) 
as representative for the loss functions at small q a b and 
arbitrary values of q c . In this way, we obtained a good 
fit to the measured self-energy in the optimally doped 
cupratea*. 

B. a 2 F(ui) with metallic screening effect 
(overdoped) 

In the following, we will include the metallic screening 
effect to calculate i?e(S) in the overdoped case. In the 
experimental paper— we found that the self-energy in the 
overdoped system was reproduced in detail, 

a 2 F(uj) = b^^- ? 7m ( ^] (12) 

using the loss function of the optimally doped system and 
a 'filter' function representing a characteristic 'screening' 
frequency Lo S rc,c — 60meV. As we will now argue, the 
'filter' function uj 2 /(ui 2 rcc + lo 2 ) is consistent with the 
changes in the screening seen in optical spectroscopy. To 
incorporate the effect of the metallic screening, one could 
model the c-axis polarizability according to Eq. (JTJ as, 

uj 2 

Vq//c P e ( q//c ^) = - T f^- f (13) 

We have neglected the details of the frequency and 
fermionic momentum dependence of the scattering rate. 
This simple Drude form is an over simplification but all 
what matters for the phonon screening are the overall 
energy scales of screening and the phonon frequencies. 
Moreover, we use "damping" as a way of loosely pa- 
rameterizing the change from coherent metallic transport 
along the c— axis in overdoped systems to incoherent in- 
sulating behavior near optimal doping. The nature of 
this change is a very important issue, relating to charge 
confinement and/or the incoherent nature of anti-nodal 
quasi-particles governing c— axis propagation in an LDA 
treatment. A full description of this change is difficult 
around optimal doping where strong correlations are de- 
veloped. An adequate description of the role of corre- 
lations requires treating both electronic correlations and 
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EP interactions on equal footing. Since such treatments 
are presently limited, wc give only a heuristic picture of 
the consequences of strong damping and incoherence. 

For finite u> PtC , l/e e ; will also acquire an imaginary part 
which contributes to the Eliashberg function. When the 
plasmon is underdamped while uj p _ c is in the phonon fre- 
quency window, the plasmon would become well 'visi- 
ble' in the self-energjiit However, in the overdamped 
regime, T > to P}C , of interest to the c-axis of the cuprates, 
this turns into a smooth background contribution. In 
the experimental determination of the self-energy such 
backgrounds are subtracted. Although it would be in- 
teresting to find out if such backgrounds associated with 
the c-axis metallization can be extracted from the ex- 
perimental data, these should be ignored in the present 
context and we should focus on the phonon contribution 
to the effective interaction in Eq. ©, 



W sc - ph (q,uj) 



M 2 q D(q^) 
e 2 (<Z» 



(14) 



where e(q, oj) = e e i/eoo and e e i is modeled according to 
Eq. Q). The renormalized phonon propagator— is, 



D(q,uj) = 



1 - Mp (5»P e ( 9 »/e( ? » 



(15) 



By using Eq. (0 and the same assumptions as for the 
optimally doped case, the Eliashberg function including 
the screening effect becomes, 



a 2 F(w) 



\M. 



9 ~ol 



ImD{u) 



(16) 



This more general form does reproduce the self- 
energies of both the overdoped and optimally-doped sys- 
tems quite well, as long as the c-axis charge dynamics is 
not turning strongly underdamped. 



Given the a 2 F(uj) as in Eq. l|16f) . the factor 



iR 1 



acts 



like a filter function in the overdamped regime. The char- 
acteristic screening frequency 0J S cr,c can be determined by 
demanding — y| 2 = 1/2 such that the phonons with 
frequencies less than uj scr can be regarded as screened. 
Since e(w) = e e z(w)/eoo it follows immediately from Eq. 
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4w„ 



1/2 



(r 2 + 2J V 



1/2 



where oj' p c = 
(T >> uj p ,c), w s 




In the overdamped regime 
, iC as in Eq. ©. 



This cj scr ,c can be directly related to the electronic 
background in the optical conductivity, cr^(w). Eq. 
leads to a simple Drude form for the conductivity which 
can be expanded in the regime r >> uj as, 



(Tl (Oj) 



4ttT [ 



r 2 



+ ...) 



(18) 
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FIG. 1: The experimental c-axis optical conductivity in op- 
timally doped (OP) from Ref. 16, (h). Various models for 
the optical conductivity of overdoped (OD) Bi2201 (e)-(g). 
The corresponding Eliashberg functions (solid line) and elec- 
tronic self-energies (dash line) assuming that the predomi- 
nant electron-phonon coupling is due to the unscreened c-axis 
phonons (a)-(d); experimental data in symbol are from Ref. 
4. Parameters: (a)-(c), u>scr,c ~ 60 meV. (a) T < lu p , c , T — 45 
meV and lj p , c = 190 meV, (b) F > u p ,c,T = 430 meV 
and LOp.c — 350 meV, and (c) T w 



p ,c,r = 180 meV and 
= 250 meV. (d) OP case, T = 70o' meV and uj p , c = 250 



meV. 



and it follows from Eq. [21 that the height of this back- 
ground corresponds to w 2 jC /47rr = (e o/47r)ix' scriC in the 
overdamped regime (r >> lo p , c )- One notices that in 
principle the self-energy could be determined using the 
information in the optical conductivity, avoiding any free 
parameters. We note that since the optical data on over- 
doped Bi2201 are not available and we model these ac- 
cording to what has been measured in the La2- a: Sr :l: Cu04 
(LSCO) system. 



III. COMPARISON WITH EXPERIMENTAL 
DATA 

According to Ref. 9, the ab plasma frequency u> Pt ab ~ 1 
eV while w Pi0 fc/ ^ P , c ~ 4 to 10 and we take w PjC ~ 
^p,ab/4 ~ 250 meV. To reconstruct the measured optical 
conductivity in the OP case 16 [Fig. 1(h)] we need a large 
r ~ 700meV to obtain a background w 12 Sl~ 1 cm~ 1 . In 
Fig. l(f)-(g) we combine the (dressed) phonon spectrum 
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of Bi2201 with the background as measured in the LSCO 
system at 30% dopingil, translating in an u> scr . c ~ 60 
meV, indicating the insensitivity to the precise values of 
bj p _ c and r separately. In Fig. 1(b)- (d) we show the 
outcomes for both the Eliashberg functions and the cal- 
culated self-energies, finding results which closely track 
the phenomenological filter function used in Ref. 4. 

For completeness we include the outcomes assuming a 
strongly underdamped plasmon [Fig. 1(a) and (e)]. One 
sees immediately that this yields a much less satisfactory 
outcome. We also notice that in the overdamped cases 
Fig. 1(b)- (d) the 'phase space' parameter b ~ 1 in all 
cases while it has to be strongly reduced assuming the 
underdamped plasmon (b ~ 0.35). 



IV. CONCLUSION 

We have demonstrated that the large scale changes 
found in the self-energy of the nodal quasiparticles in 
Bi2201 are in detailed quantitative agreement with the 
measured changes in the screening properties along the 
c-axis when it is assumed that this self-energy is mostly 
due to the scattering of the electrons against the polar 
phonons associated with the motions of ions in the c- 
direction. This is a direct evidence for the presence and 
importance of this type of interaction and in a future 
publication we will elaborate further consequences of this 
unconventional electron-phonon coupling 15 . 



APPENDIX A: EFFECTIVE PLASMON 
FREQUENCY, w tA 

In a coupled layered electron gas, the dielectric func- 
tion using the RPA approximation (Eq. (7) in Ref. 9) is 
given by 



APPENDIX B: EFFECTIVE INTERACTION, W, 



S 



e(w,£?) = 1 - 

q 2 ' \ u,' ^' 



lu(lo + i0 + ) q 
1c 2 o f v cQc v F q ab \ 
a V u) w J 



where S is a form factor, v c is the effective velocity per- 
pendicular to the layers, and p° is the function 



p°( a ,b) = — [ {(l-acos0) 2 -b 2 \ 
™ Jo 



1/2 



cos 



And, hence the effective 3D plasmon frequency oj v ^ 
may be written as 



off 



The detailed derivations of the following can be read 
from Ref. 13. When considering the single-process scat- 
tering of two electrons: i) by the unscreened electron- 
electron interaction and ii) by sending a phonon from 
one to the other, the combined interaction W° can be 
written as: 

W°(q,uj) = v™ + V ph (q,uj) 

where = Aire 2 / (q 2 e OQ ) is a bare Coulomb part and 
V p h is a phonon mediated part. 

This phonon mediated part (same as Eq. (|SJ ) is given 
by, 



V ph (q-,Lu) = M 2 D"(q,u) 

where M q is a bare EP vertex and D° is a bare phonon 
propagator. 

Next, when considering all the possible of multiple- 
process scattering, the effective interaction in the classic 
RPA form, W e //, can be written as: 

W°{q,u) v q 
effKQ, ) 1 _ W o^ u))Pe ^ )U) e tot %u>) 

where e tot (q,oj) = - v q P e (q, u) + e ph (q, u), P e (q,uj) is 
an electronic propagator and e p h(q,uj) is a phonon polar- 
izability given by Eq. Q). 

To explicitly separate the electron-electron part of the 
effective interaction, W e ff can be rewritten to be the 
same as Eq. ( or: 



,, 2 q < 
q 



lb,, 2 n (VFqab v c q c \ q\ 2 2 Q fv c q c 

T^p.abP [ 1 ) "I 2 W P,cP I ! 



VFQab 



OJ OJ 



e tot (q,oj) e ei (q,u;) 



sc—ph 



where W sc - P h is the screened EP interaction given by Eq. 
(|14|l and the renormalized phonon propagator is given by 
Eq. fT5). 



APPENDIX C: NUMERICAL FACTOR, b 

In the insulating case (e e ; = too), from Eq. Q and 
(| 1 1 P . b can be written as 

l r d 2 k' 



(2tt)3 



I^r\M n _ g ,\ 2 ImD°(k-k',oj) 



I M 2 

l -^^ImD°(oj) 

In the overdoped case, replace M with e^M/e^ and D° 
with D in the above equation. 

When momentum dependences are weak, the factor b 
can be approximated by 



where in the limit q — > 0, S — > 1 and p — > 1. 



(2tt) 3 



Vp 
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